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Two issues are treated in this work: (i) the generic fact that if a fermionic superfluid in the BCS 
regime overflows from a narrow container into a much wider one, pairing is much suppressed at 
the overflow point. Physical examples where this feature may play an important role are discussed, 
(ii) A Thomas-Fermi (TF) approach to inhomogeneous superfluid Fermi-systems is presented and 
shown that it works well in cases where the Local Density Approximation (LDA) breaks down. 



Superfluid fermions in finite systems can exist in traps 
of cold atoms, in nuclear systems, in small metallic clus- 
ters, etc. An interesting question arises what happens 
to the superfluid if its Fermi level reaches the edge of a 
finite container, i.e. either the fluid overflows into the 
continuum or it pours into another container of much 
larger dimension. A trapping potential of this type has 
experimentally already been generated for the study of 
cold bosonic atoms It should also be possible to use 
it for fermionic atoms (ij. 

In the inner crust of neutron stars, there also may occur 
the situation where coexists a superfluid neutron gas of 
variable density in between the lattice of (superfluid) nu- 
clei [3, 15 . This situation often is mimicked by a Wigner 
Seitz cell, where the single particle potential has a pocket, 
representing the nucleus, embedded in a large container. 
Still other systems may exist with similar situations. 
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FIG. 1: Schematic view of the potentials used in this work. 
Panel (a) shows a prespective view of the slab potential which 
is translationally invariant in x, y direction. Panel (b) repre- 
sents a spherically symmetric container 



The purpose of this work is to study superfluidity 
of fermions at the overflow (drip) in the BCS regime. 
Since the quantal solution of BCS equations in geome- 
tries with rapidly varying single particle potentials with a 
large number of particles is numerically difficult, we will 
present, as a second objective of this work, a Thomas- 
Fermi approach to inhomogeneous superfluidity which 
shows good performance in situations where LDA fails. 

For our study we first will use a schematic model of 
slab geometry with a transverse potential of large ex- 
tension L posessing at the origin a 'pocket' of variable 
depth and size R much smaller than the outer container. 
Schematically such a potential is shown in Fig. la. This 
slab configuration may roughly mock up one sheet of a 
so-called Lasagne configuration in the inner crust of neu- 
tron stars fHI]. We, therefore, will use nuclear dimensions 
for the slab model but they can easily be replaced by di- 
mensions relevant for other systems. Our model and the 
ensuing generalisations treated below, therefore, are be- 
lieved to be generic. We will study the slab configuration 
also because the quantal solution of the gap equation is 
evaluated relatively directly and the quality of the TF 
approach can thus be established. Once this is achieved, 
we also will go over to other geometrical configurations. 
We, for instance, will treat a second potential shown in 
Fig. lb with spherical symmetry, a kind of which, as al- 
ready mentioned, has been used for bosonic atoms in pO] . 

The wave functions and eigenenergies of a box as shown 
in Fig. la with a potential-hole are given in For pair- 
ing, we use a contact force with a cut off A, to make things 
simple. Integrating over momenta in slab direction the 
usual gap equation A„ = Y.n' I (■Lny^ '^nn'Qi^ - £n' - 

£p)A„V(2£:„'(p)) with Enip) = v^(£„ + £p - At)2 -t- A2 

the quasiparticle energy, Q{x) the step function, and 
Sm£p being the discrete single particle energies in trans- 
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FIG. 2: (Coloronline) Quantal and TF pairing gap in the slab 
geometry as a function of the chemical potential. 



verse direction and kinetic energies in slab direction, re- 
spectively, one arrives at the following gap equation 



with 



(1) 



FIG. 3: (Coloronline) Position dependence of the gap in the 
slab geometry for different values of the chemical potential. 
Quantal, TF, and LDA results are shown. Notice that A for 
= 0.5 and -5.0 MeV is practically zero in the gas region. 



g (E) is the corresponding level density to lowest order 
in /i, i.e. the usual TF expression Quantum numbers 
and energies are simply characterised by the continuous 
energy variable E which takes the place of the discrete 
values £„ in the quantal case. 

The TF version of the gap equation (1) then reads 



= ^ A„ In (2) 



where m is the particle mass and the indices n stand 
for the level quantum numbers in the confining po- 
tential of Fig. la. The matrix elements Vnn' = 
—g J^l^ \ipn{z)\'^\(pn'{z)\'^dz of the contact force —gS{r — 
r') can be evaluated straightforwardly from the wave 
functions (/?„(z) given in Q. 

For an example we take as cut off A = 50 MeV counted 
from the edge of the pocket potential whose depth be 
Vq = —40 MeV. Its extension ranges from —R to +R 
with i? = 10 fm. The wide potential with infinite walls 
has extension from —L to +L with L = 100 fm. For the 
coupling strength we take g= 150 MeV fm'^. 

Before we show the results, let us explain our Thomas- 
Fermi (TF) approach for this problem. For this, we 
transform Eq.(l) into a continuum version in the fol- 
lowing way. We first consider the Wigner transform 
of the density matrix corresponding to the state \n): 
[Pn]w = [|'^)(lillw/ and take the fi, — >■ limit of this ex- 
pression, see 



fE{z,p) 



1 



5{E^H,i) + 0{h') 



(3) 



where Hd. = ^ + ^(^) is the classical Hamiltonian with 
V{z) the potential (Fig. la) in transverse direction and 



A{E)^- [ dE'g{E')V{E,E')K{E') 

Jvo 



(4) 



with K{E) an obvious generalisation of Kn- The matrix 
elements V{E, E') can be evaluated in replacing |(/3„(z)|^ 

by ^ pW{z) = -^^J-{^^f'[E-V{z)]-^/^ the on- 
shell TF density in transverse direction. As the reader 
will easily realise, the way of proceeding is very different 
from usual LDA where the finite size dependence is put 
into the (local) chemical potential whereas here it is put 
into the matrix elements of the pairing force. 

We are now in a position to solve the quantal and TF 
gap equations for the above mentioned parameter values 
of our model. The result for the gap at the chemical po- 
tential /i is shown in Fig. 2 as a function of fi. We start 
with fi from the bottom of the pocket well, i.e. with zero 
density. We then increase fj,, i.e. the density. We see that 
once the fill up of the pocket reaches its top, the values of 
the gap sharply drop and practically reach zero. In the 
continuum the gaps slowly rise again. We see that quan- 
tal and TF values are in close agreement. The overshoot 
of the TF solution for negative ^ very likely is due to 
the smallness of the pocket which only can accomodate 
nine bound levels. It may be partially cured including h 
corrections Q which, however, we do not consider here. 
Before we come to an explanation of the drop of the gaps 
at overflow (drip), let us study the gaps as a function 
of position in transverse direction: A(z) = —gK{z) with 
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Semiclassically, this expression 
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FIG. 4: (Coloronline) Average TF gaps at the Fermi energy 
as a function of the chemical potential for the potential shown 
in Fig. lb. In the completely filled optical trap (/x = V) we 
accomodate 10^ atoms in each spin state. The total number of 
atoms in the trap with ^/huiopt='^Q, 80 and 120 are indicated 
in the upper horizontal axis. 
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FIG. 5: (Coloronline) Upper panel: Radial dependence of the 
TF gap in the considered WS cells. The end points indicate 
the radius of the WS cells. Lower panel: Comparison between 
TF and LDA gaps as a function of the position in a WS cell 
containing a single to^Zr nucleus. 



tdEg^^{E)K{E)pY{z). 



K{z) = Y^KnWn 

becomes: K{z) — jy^' 

In Fig. 3, we show the density profiles for three values 
of /i: jl — 40, 0.5, and - 5 MeV. We see that quantal 
and TF results agree, up to shell fluctuations, very well. 
We also show the LDA results. We see that they can be 
as wrong as by 50 percent. For other choices of system 
parameters the LDA error may even be worse. This stems 
from the fact that in TF (and, of course, also quantally), 
there is coupling between inside and outside the pocket, 
i.e. the Cooper pair wave function extends into both 
regions what tends to equilibrate the values of the gaps. 
In LDA the contrast is much too strong. The drop of the 
gaps when crossing the threshold can be explained by the 
fact that the single particle states are strongly delocalised 
in the outer container and, thus, their contribution to the 
pairing matrix element Vn^n' becomes very small. 

Having gained faith into our TF approach, we now 
can explore other geometries and other systems, more 
difficult for quantal solutions. In Fig. 4 we show the 
result for A in the spherical double harmonic oscilla- 
tor potential shown in Fig. lb which may be realised 
with cold fermionic atoms to study the overflow situa- 
tion. A zero range force with strength g=-l.Q and cut off 
A=164.34 (in the corresponding optical trap units with 
ujopt = 27r X 1000 Hz taken from [1]) is used. We see that 
the result is qualitatively similar to the slab case though 
in this spherical geometry the dip does not quite reach 
zero and also is shifted slightly to an energy above the 
break of the potential. Note that this depends strongly 
on the choice of the ratio (jJmag/^opt as it can be seen 
in the figure. Also the gap starts to decrease towards 
the minimum quite early. It shall be interesting to see 
whether our prediction can be verified experimentally. 

Let us now make a more realistic study of Wigner-Seitz 
(WS) cells including electrons in /3 equilibrium to simu- 
late the inner crust of neutron stars [5|. To this end the 
mean-field is computed selfconsistently using the BCP 
energy density functional [l^l together with the TF ap- 
proach as explained in 11| . The semiclassical description 
of the WS cells including pairing correlations at TF level 
is obtained from this mean-field using the finite range 
par t of the Gogny DIS force [l3l in the pairing channel 
[13l |. It must be pointed out that the total energy per 
baryon obtained with our TF approach is in very good 
agreement with the old quantal calculation of Negele and 
Vautherin IJ] as it is explicitly discussed in Ref. ll|. In 
Fig. 5 are displayed the corresponding gaps with their 
radius dependences. It is seen that when the gap is small 
outside the region of the nucleus, then the gap also is 
small inside the nucleus. This stems from the very large 
coherence length where one neutron of a Cooper pair can 
be in the huge volume of the gas and the other inside the 
small volume of the nucleus (proximity effect). In this 
way the gas imprints its behavior for the gap also inside 
the nucleus. Such a conclusion was also given in a quantal 
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Hartree-Fock-Bogoliubov (HFB) calculation by Grasso et 
al. in [isj what shows that the here employed BCS ap- 
proximation apparently yields very similar answers as a 
fuU HFB calculation for WS cells BEI- Finally, in the 
lower panel of Fig. 5 we show a comparison of LDA and 
present TF results for a particular WS cell. We see a 
huge difference in the surface region of the nucleus. This 
simply stems from the fact that in this case of the 49° Zr 
nucleus in the WS cell the gap is very small and, there- 
fore, the coherence length very large invalidating LDA. 
A study with examples a little less unfavorable for LDA 
is given in [17]. 

For isolated nuclei at the neutron drip the situation 
may be different. It seems that in this situation the dif- 
ference between HFB and BCS approaches may be signif- 
icant. Somewhat conflicting results in this respect exist 
in the literature. In ref [l8| very similar results to ours 
are found for S-wave pairing. On the other hand in [Toj 
the gap seems to rise towards the drip before it bends 
down. Similar results have recently been found in [20| . 
Preliminary investigations show that these discrepancies 
may be due to large shell fluctuations in isolated nuclei. 
More studies in this direction seem to be necessary. 

Summarizing, we have studied superfluid fermions in 
a large container, either external (cold atoms) or created 
self consistently (nuclei) for situations where the top of 
the fluid reaches the edge of a small pocket situated at 
the origin of the wide confining potential. The gap drops 
to zero at the edge before rising again when the density 
fills up the outer container. This at first somewhat sur- 
prising phenomenon can be explained quite straightfor- 
wardly. Such situations, as already mentioned, can exist 
in cold atoms and nuclei in the inner crust of neutron 
stars, two examples treated here with their specific form 
of containers. For small systems like isolated nuclei at 
the neutron drip, the situation may be blurred by shell 
effects. 

As an important second aspect of this work, we showed 
that a novel Thomas-Fermi approach to inhomogeneous 
situations can cope with situations where LDA fails. This 
means that our TF approach is free of the restrictive con- 
dition, prevailing for LDA, that the Cooper pair coher- 
ence length must be shorter than a typical length I (the 
oscillator length in the case of a harmonic container) over 
which the mean field varies appreciably. On the contrary, 
our TF theory has the usual TF validity criterion, namely 
that local wavelengths must be shorter than I. 

The accuracy of our TF approach opens wide perspec- 
tives for a treatment of inhomogeneous superfluid Fermi- 
systems with a great number of particles not accessible 
for a quantal solution of the BCS (HFB) equations. Such 
systems may be cold atoms in deformed containers (even- 
tually reaching millions of particles), superfiuid-normal 
fiuid (SN) interfaces, vortex profiles, etc. As a matter of 
fact, as is well known the TF approach becomes the 
more accurate, the larger the system. Thus the TF ap- 



proximation is complementary to the quantal one in the 
sense that the former works where the latter is difficult 
or even impossible to be obtained numerically. 
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